The transition between hole-pairs and four-hole clusters in four-leg t-J ladders 



(N 
O 

o 

(N 
(N 



I 



I 

o 
o 



> 

cn 
o 

O 



I 

o 
o 



X 
J3 



Thomas SiUer \ Matthias Troyer\ T. M. Rice\ StevenR. White^ 

1 Theoretische Physik, Eidgendssische Technische Hochschule, CH-8093 Zurich, Switzerland 
^ Department of Physics and Astronomy, University of California Irvine, CA 92697 

(February 1, 2008) 

Holes weakly doped into a four-leg t-J ladder bind in pairs. At dopings exceeding a critical doping 
of 5c — ^ four hole clusters are observed to form in DMRG calculations. The symmetry of the 
ground state wavefunction does not change and we are able to reproduce this behavior qualitatively 
with an effective bosonic model in which the four-leg ladder is represented as two coupled two-leg 
ladders and hole-pairs are mapped on hard core bosons moving along and between these ladders. At 
lower dopings, 5 < 5c, a, one dimensional bosonic representation for hole-pairs works and allows us 
to calculate accurately the Luttinger liquid parameter Kp, which takes the universal value Kp — 1 
as half-filling is approached. 



I. INTRODUCTION 

Strongly correlated electrons confined to move on lad- 
ders formed from coupled chains have been an active 
topic of research in recent years. These examples of 
lightly doped spin liquids can be very efficiently ana- 
lyzed using the density matrix renormalization group 
(DMRG) method, which allows long ladders to be 
accurately simulated. Ladders with an even number 
of legs have spin liquid groundstates at half-filling and 
the spin liquid character remains upon doping due to 
the binding of holes to form pairs at low concentrations 
P . The underlying physics of the ladder systems has 
close similarities to that of a doped resonant valence bond 
(RVB) phase which is widely believed to be realized in 
the high-Tc superconductors. The numerically accurate 
DMRG simulations of the strongly correlated t-J model 
PI give a clear picture of the behavior of hole-pairs but 
only in qualitative terms as the density varies. 

In this paper our aim is to analyze the DMRG re- 
sults and to use them to determine the effective inter- 
actions between hole-pairs which govern their behavior. 
Since the spin sector of even leg ladders remains gapped, 
there are only charge degrees of freedom at low energies 
and these can be represented in terms of suitable hard 
core boson models. These in turn can be analyzed much 
more easily and completely. In an earlier paper jl^ we 
described such a representation or mapping of the hole- 
pairs on a two-leg ladder to hard core bosons on a single 
chain. This enabled us to determine the effective repul- 
sive interactions between the hard core bosons by fitting 
to the hole density distribution determined by DMRG 
for the t-J model. In the present paper we extend this 
analysis to four-leg ladders. The extra transverse de- 
grees of freedom allow the formation of larger clusters 
of two hole-pairs (i. e. four hole clusters) when the hole 
density exceeds a critical concentration 6c — ^- A clus- 
ter represents not a bound state of the hole-pairs but a 
finite energy resonance which becomes populated when 
the chemical potential reaches a certain threshold. 

The outline of the paper is as follows. First we treat 
the low concentration region where the density profiles 



obtained from DMRG simulations show that hole-pairs 
simply repel each other. In this region we extend the 
mapping onto the model of hard core bosons that we 
used earlier for the two-leg ladder and obtain a parame- 
terization of the repulsive interactions and Luttinger liq- 
uid exponent Kp. We compare the evolution of Kp with 
hole density in the two and four leg ladders. In both cases 
Kp ^ 1 as the hole concentration 5 — > indicating pre- 
dominantly superconducting correlations. To treat the 
formation of the four hole clusters it is necessary to in- 
troduce additional transverse degrees of freedom. This is 
done by mapping the hole-pairs on a four-leg ladder onto 
a hard core boson model on a two-leg ladder. Choosing 
a potential with a repulsive tail but with a finite energy 
resonance on a single rung allows us to reproduce the 
DMRG results for the four-leg t-J model, both for the 
kink in the chemical potential and for the change in the 
density profile at the critical density. The paper ends 
with a concluding section. 



II. HOLE PAIRS AND CLUSTERS 

We consider the t-J ladder Hamiltonian 

H = -t '^(^.Cla + cl,C,,.)V (1) 



+J ^(SiSj - jniUj) 



<i,j> 

for a four-leg ladder (4LL) configuration as depicted in 
Fig. |l]. Here cr — (f , j) denotes the spin index and (i, j) 
the summation over nearest-neighbor sites with i = (x, y) 
as site index. The operators and Ci^„ create or de- 
stroy electrons with spin a at position i respectively. The 
corresponding density operator is ni^a ~ c} cr^i,<j and Ui — 
rij I -I- ni^|. The projection operator V = Yiii^ ~ "-i,?"^!.!) 
prohibits double occupancy of a site and Si denotes the 
spin at site i. We use J = 0.35 1 throughout. If not noted 
otherwise, we use the DMRG method to compute ground 
state properties which implies that open boundary con- 
ditions on the ends of the ladder are used. 
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FIG. 1. Indexing convention for a (L x 4) ladder. 



We start by examining the internal structure of a hole- 
pair (HP) on a t-J 4LL, by studying the hole-hole corre- 
lation (HHC) function. Since the DMRG computations 
are performed with open boundaries, special care must 
be taken in the measurement of the HHC-function, which 
we define as 



9y,y 



{<,v){<',y) 



(2) 



Here vJ^ y = 1 — nx,y denotes the hole density operator 
acting on position (x, y) and n'' is the average hole den- 
sity. Introducing relative and center of mass coordinates, 
r = (0, 1, . . .) and Rq = (1, 3/2, 2, . . .) in the long direc- 
tion, we define with 



(3) 




FIG. 2. Hole-hole correlations g{r) as a function of the 
relative separation r along the legs measured on a (32 x 4) 
ladder with two holes and J = 0.35 1. The inset shows a 
logarithmic plot of the same data. 



the HHC between different rungs. To obtain a good ap- 
proximation for the HHC function g{r) on an infinite 
ladder, we measure g{r, Rq) for an {L x 4) ladder in the 
middle of the system, i. e. for Rq ~ L/2. In this way we 
have measured g{r) on a (32x4)-ladder with open bound- 
aries. The result is plotted in Fig. ||. It shows clearly that 
in the ground state the two holes are bound. The cor- 
relation function decreases exponentially as can be seen 
from the logarithmic plot in the inset. The correlation 
length is 1.358 and hence larger than in the two-leg lad- 
der (2LL) case where we have obtained 1.184 jl^. This 
suggests that the binding energy of a HP is reduced as 
the effective dimension is enhanced. With E{N) as the 
ground state energy with N holes doped into the ladder 
we have computed the binding energy defined as 



Eb = 2E{1) - E{0) - E{2) 



(4) 



for two and 4LL using systems with 40 and 24 rungs re- 
spectively. As expected, we have obtained a lower bind- 
ing energy Eb = 0.1062 1 for the 4LL than for the 2LL 
case, Eb = 0.1487 i. 

If the HP behave as hard core bosons (HCB) at low 
hole doping one should find as many maxima in the hole 
density profile of a t-J 4LL with open boundaries as the 
number of HP. In fact, we observe this only below a crit- 
ical doping of Sc — |_in contrast to the results obtained 
for two-leg ladders, |l2| where this equality holds to much 
higher dopings. Figures || and ^ shows the hole density 
profiles for the groundstate of a (16 x 4) ladder for 2 up 
to 12 holes. The plot with 6 holes doped into the ladder 
corresponds to a hole density per site <5 < | and shows 
three maxima as expected. For 8 holes doped into the 
ladder (5 = i) we could expect to find four maxima, but 
we observe only three, suggesting that larger clusters of 
four holes are forming in some way as will be discussed 
later. 
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FIG. 3. Hole density profiles {n'^,y) on the inner and outer 
legs of a (16 x 4) t-J ladders with J = 0.35 1 for Nh = 2, 4 
and 6 holes doped into the ladder (5 < |). 



FIG. 4. Hole density profiles (n^^y) on the inner and outer 

legs of a (16 x 4) t-J ladders with J = 0.35 1 for Nu = 8, 10 
and 12 holes doped into the ladder {5 > i). 
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FIG. 5. Mapping of hole-pairs from a four-leg ladder to 
hard core bosons on a chain. The center of mass coordinate 
R of the hole-pair determines the hard core boson position. 

III. LOW DENSITY OF HOLE-PAIRS 

In this section we examine the t-J 4LL at low hole 
doping S, i.e. S < Sc- In order to set up an effective 
bosonic model we first check the symmetry of the ground 
state at low hole dopings using a Lanczos method to get 
the lowest lying energy states with exact diagonalisation 
for a 6 X 4 t-J ladder with periodic boundary conditions 
along the legs. For 2 and 4 holes doped into the half filled 
ladder the ground state has even parity under reflection 
along and perpendicular to the legs. 

Since the spin part of the ground state wavefunction 
is a singlet and thus antisymmetric and the spin excita- 
tions are gapped, Q the charge degrees of freedom can 
be described by an even parity wavefunction where HP 
are considered as effective (hard core) bosons. At low 
hole dopings, S < Sc, only one HP is found on a rung. 
Hence, we can map HP in the {L x 4) t-J ladder with pe- 
riodic boundary conditions in the long direction on HCB 
on a single closed chain as shown in Fig. ||. This effec- 
tive model is completely analogous to the model for 2LL 
introduced in Ref. 

Since holes on a ladder can pair with more weight on 
adjacent rungs, our effective model incorporates the pos- 
sibility that the "center of mass" of a HP can lie on a rung 
or between two rungs as shown in Fig. ^ Note that for 
even (odd) distance r along the legs between two holes, 
the center of mass lies on a rung (between two rungs). 
The HHC function g{r) is connected to the probability 
of finding the center of mass of a pair on a rung, Wint , or 
between two rungs, Whaif , by 



Wint 



(5) 



I-odd 



The same occupation probabilities can be obtained 
with the one boson Hamiltonian 



Hb 
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-t* E (B^B, 



-Bl^M- 
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Nr (6) 



for a boson which moves on a closed chain with length 
L, i. e. 2 L sites, under the action of an alternating on-site 
potential which is or e on integer or half integer sites 
respectively. Here Nr = bI^Br and R= (1/2, 1, ... , L). 

Figure || shows the mapping of HP from the ladder to 
effective bosons on a single chain. Note that the center 
of mass coordinate R of the pair determines the position 
of the effective boson. Here the ratio of the probabilities 
If half /wint to find the boson on a site with half integer or 
integer R depends only on e/t* and can easily be obtained 
as 



Whalf 



+ (4<*)2 -e^^TTiFp 



Wint £2 + {4t*y + eVe^ + (4F)2 
From this equation and Eq. (^ we obtain for e 



(7) 



4 1* sinh 
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I^rodd air) 



(8) 



With g{r) = g{r, Rq) obtained as explained in Sec. || for 
i?o — L/2 from a (32 x 4) t-J ladder we get Whaif/wint — 
1.18 and e = -0.33 i*. This result shows that the HP is 
mainly centered between two rungs. 

Once the boson density (Nr) for the model (^) has 
been computed, we obtain the hole rung density (nj) — 



J2y = l 



p!^ y) on rung x by the convolution 



i S.„,,g(0((^x-j) + (A^.+r)) 

Next we generalize the one boson model to a finite 
density including the interactions between the HCB and 
also the effect of open boundary conditions and write the 
effective Hamiltonian for the HCB as 



(10) 



The potential Vi^t gives the interaction between HCB, 
i. c. HP in the t-J model. Since the simulations are for 
finite systems, we have to take into account the interac- 
tion of the HP with the boundaries. The potential Vh 
has been introduced to describe this effect. 



A. Computing Vh and Vint 

Our procedure is to compute the HCB density (Nr), 
convolute it with g{r) according to Eq. (||) and then to 
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FIG. 6. Hole rung density (n^) for two holes on a (24 x 4) 
t-J ladder with J = 0.35 1 computed directly and with the 
effective model using smooth boundary conditions. For the 
effective model the data for = and the best fit is shown. 

compare it with the hole rung density (n^) of the corre- 
sponding t-J system. In this way we obtain Vb and T^nt 
by fitting the density profile of the effective model to that 
of the t-J model. 

However there is a problem in the density profiles of 
the t-J ladders. The open boundary conditions induce 
density oscillations at the boundaries which cannot be 
neglected as in the case of the 2LL. Due to the 
higher number of degrees of freedom the HP in the 4LL 
gets more distorted, when it approaches the open ends of 
the system. To circumvent this problem we use smooth 
bounday conditions (SBC) as proposed in Ref. for 
both, the t-J and the effective model. The SBC introduce 
smoothly decreasing energy parameters into the Hamil- 
tonian as the open ends of the ladder are approached. We 
use SBC only in the long direction which extend into the 
t-J ladder until the fifth rung from the open ends. We 
choose the same smoothing function as used in Ref. 
for the Hubbard model. Details are explained in Ap- 
pendix 

We obtain Vb by considering one HP in the t-J model 
and choose an exponentially decreasing form for the po- 
tential term Vb 

Vb = Vb2_^NR [e 56 -He j (11) 

R 

with the two parameters Vh and ^b- Figure ^ shows the 
results of the fit and the optimal parameter values are 
displayed in Table I. Note that the rather large value of 
Vh is a consequence of the smoothing function, which also 
decreases the energy parameters of Vb, i.e. Wb, near the 
boundaries. 

To obtain the interaction potential T^nt we proceed in 
the same way as for Vb and choose a hard core form 

Vnt = 51 II «int(l^ - R'\) NrNr. (12) 
R R'>R 



TABLE I. Parameters for the boundary potential Vb given 
in Eq. (^l|) obtained from the density profile of a (24 x 4) t-J 
ladder for J = 0.35 1. 
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FIG. 7. Hole rung density (n^) for J = 0.35 1 calculated 
for four holes on a (24 x 4) t-J ladder computed directly and 
with the effective model using smooth boundary conditions. 
For the effective model the data for v^nt = and the best fit 
is shown. 

{OO Ai? < Aiimin 

vi AR = AR^in . 

„g-(i?-i?„.„-i)/« Ai?>Ai?,ni„ 

We consider a (24 x 4) t-J ladder with four holes and 
the corresponding effective model with two HCB and use 
three fit parameters for the interaction potential and the 
additional parameter ARmin, which broadens the hard 
core of the bosons, to model the interaction between the 
pairs. 

Using the fits to the density profiles as shown in 
Fig. 1^ the parameter values quoted in Table |l| are ob- 
tained. Figure § shows the result together with that from 
Ref. [Hi for the two leg ladder. 

We have tested the results by comparing the density 
profiles for various numbers of HP and ladder lengths. 
We find good agreement between the density profiles ob- 
tained from the t-J model and the effective model for 
hole dopings S < 0.08, as can be seen from Fig. ^. 

B. Luttinger liquid parameter Kp 

Having determined the effective HCB model that de- 
scribes the low energy properties of the HP in the t-J 
model, we can express the charge density and the super- 
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FIG. 8. Hole rung density (n^) for J — 0.35 f for six holes 
on a (32 x 4) t-J ladder using smooth boundary conditions 
and computed with the interaction potential obtained from 
fits with two hole-pairs on 24 x 4 sites. The data for Vint ~ 
for the effective model are also shown. 



TABLE 11. Parameters for the interaction potential Vint 
given in Eq. (^) obtained from the density profile of a (24 x 4) 
t-J ladder for J = 0.35 1. 
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FIG. 9. Interaction parameter Vint/t from Eq. |l^ obtained 
for the four-leg ladder in comparison with the result from the 
two- leg ladder from Ref. [ p^ . 



i4 



0.6 



0.4 



- four leg ladder 
■a- two leg ladder 




0.05 



0.1 



0.15 



0.2 



0.25 



FIG. 10. Luttinger liquid parameter Kp as a function of 
hole doping for the four-leg t-J ladder in comparison with the 
results for the two- leg ladder from Ref. The hole doping 
is given in holes per rung, 5r = Ni S, with A*'; as the number 
of legs. 



conducting correlation functions in terms of Luttinger 
liquid parameter Kp as 



{NrNo) ~ 



const X i? ^ + const x cos{n{Nb/L) R) R 
const X i?-V2A'p 



with Nb as the number of HCB on the chain and N = 
Nb/IL the mean HCB density per site. These relations 
show that the superconducting correlations (BIBq) are 
dominant if Kp > i . For HCB in one dimension, Kp can 
be obtained from the relations mTl 



K„ 



(15) 



Here Eq denotes the ground state energy for a closed ring 
of length L, i.e. N = 2L sites, with iVb HCB and £'o(^') 
is the ground state energy of the system penetrated by 
a magnetic flux $ which modifies the hopping by the 
usual Peierls phase factor, t* i— s- t* exp(± i$/7V). From 
these two equations the charge velocity Vc can be elim- 
inated. We used exact diagonalization for HCB chains 
with lengths between 32 and 220 and with iVb ~ 2. 

The Luttinger liquid parameter Kp for the interaction 
potential given by the parameters in Table || approaches 
the universal value Kp = 1 at half filling. The results are 
shown in Fig. |l^ where Kp is plotted as a function of the 
hole doping per rung in the corresponding t-J ladder, 
5,. ^ 4:5 = 2Nb/N. For N^/N 0, corresponding to 



a very dilute HCB 
consistent with Ref. 



we have K^ 



1 



_ . Up to ,5 ~ 0.06 {5r 

superconducting correlations are dominant, 
1 

2' 



f OiN^/N) 
~ 0.23) the 
since Kn > 
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h=4, 0.05 -0.05 

FIG. 11. Real space off diagonal singlet pairing expecta- 
tion values s^j between nearest neighbor sites on a (12 x 4) 
t-J ladder for h = A. The thickness of the lines gives the mag- 
nitude. Black and shaded lines denote positive and negative 
values respectively. 

We have tested numerically the influence of the various 
parameters in Eq. ^ on Kp and Vc- A longer ranged 
interaction leads to a lower Kp and a larger Vc leaving 
the product Kp Vc nearly unchanged. As Fig. ^ shows, 
the interaction between HP on 4LL is longer ranged than 
on 2LL. Accordingly, the Luttinger liquid parameter for 
the 4LL is lower in the investigated density range. As 
can be seen from Figure |l^ the long range charge density 
wave correlations tend to overcome the superconducting 
correlations already at lower hole dopings with increasing 
dimensionality, i. e. when going from 2LL to 4LL. 

IV. THE FORMATION OF FOUR HOLE 
CLUSTERS 

In this section we examine the four-hole clusters (FHC) 
which form on a t-J 4LL at hole dopings higher than a 
critical doping 6c — ^. 

First we clarify the question whether a change in the re- 
flection symmetry occurs for the ground state wave func- 
tion with increasing doping. We consider the real space 
off-diagonal singlet pairing expectation value (SPE) de- 
fined as 

sti-{h-2\i,i,-i,i,\h) (16) 

for nearest neighbor sites Here \h) denotes the 

ground state for h holes doped into the ladder. We have 
calculated the SPE on a 12 x 4 t-J ladder for h up to 
12, i. e. a hole doping S — ■j. Figure |ll| shows the results 
for h — A. The SPE show a characteristic sign-pattern 
which would change with the reflection symmetry over 
the doping, i. e. when going from h ^ h -\- 2. We ob- 
serve no change as h is further increased and conclude 
that the ground state symmetry remains unchanged in 
the investigated density range. This result suggests that 
HP enter the same band at dopings below and above Sc- 
Figure ^ reveals that at low hole doping HP enter in 
states where they reside predominantly on the outer 2LL 
since the magnitudes of the SPE between next neighbors 
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FIG. 12. Mapping of pairs of holes from a four-leg t— J 
ladder to hard core bosons on a two-leg ladder. In this repre- 
sentation a hole-pair is centered on a 2 x 2 plaquette on one 
of the two coupled two- leg ladders. 

on these ladders are largest. But this docs not imply that 
two holes sit on the same rung. Instead the HHC func- 
tion gy yi(x,x') shows that two holes can most probably 
be found in a diagonal configuration on a 2 x 2 plaquette 
on one of the two coupled 2LL. 

To set up an effective model which can describe the for- 
mation of the FHC, we again model the charge degrees of 
freedom of HP by effective (hard core) bosons moving in 
a spin liquid but introduce additional transverse degrees 
of freedom by considering the t-J 4LL as two coupled 
2LL. This allows us to map HP - represented by their 
center of mass - from the {L x 4) t-J ladder on HCB 
on a ((L — 1) x 2) ladder as depicted in Fig. ^ In this 
representation a HP is centered on a 2 x 2 plaquette on 
one of the two coupled 2LL. The Hamiltonian is given by 

L-2 2 

H = -tiY. E(^L^»+i^^ + (17) 
1=1 j=i 

L-l 

-t,^(4li?,,2+Sl,2SM) 

1=1 

+ 1 E v{\^-^\,\J-J'\)N.,JN,,,,, 

where Bj j and Bij create and destroy a HCB on site 
respectively and Nij is the corresponding density 
operator. The matrix elements ti and tr stand for nearest 
neighbor hopping along the legs and rungs of the 2LL. 
The last term describes the interaction between HP and 
is determined by the parameter v{rx, Vy) with rx and ry 
as the relative separation of two HCB along and perpen- 
dicular to the legs respectively. 

As hopping matrix element along the legs, i/, we use 
the value obtained in Ref. (l|| for 2LL, ti = 0.303 i. The 
energy difference of the bonding and antibonding band is 
given by 2tr. We have computed this energy difference 
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TABLE III. Interaction parameter values v{rx,ry) in units 
of the hopping matrix element t of the t-J Hamiltonian (n]). 








1 


2 


3 


> 4 


vir.,0)/t 




100 


5 


0.1 





v{rx,l)/t 





0.1 


0.5 


0.1 






at wavevector — with exact diagonalisation for two 
holes on a 6 X 4 leg t-J ladder using periodic boundary 
conditions in the long direction and obtained U ~ 0.092 1 
so that ti/tr — 3. 

We are able to reproduce qualitatively the formation of 
FHC but the formation depends crucially on the form of 
the interaction between the HP. We know from Ref. 
that on a 2LL HP repel each other and hence choose a 
repulsive interaction for > which tends to zero when 
rx increases. Further we expect this repulsion to be less 
for hole pairs on different legs of the effective 2LL, i. e. 
v{rx:,l) < vlr^^O). At short distance we can expect a 
certain energy lowering due to a lower cost in magnetic 
energy in the configuration with four holes on neighbor- 
ing sites. Thus we set the constraint that v{rx, 1) has a 
local minimum at — 0. 

The solution of the two particle problem for the ef- 
fective HCB model with this type of interaction shows 
that in the ground state the two HCB are unbound. But 
there are also resonant cluster states which can be occu- 
pied at higher energies. Hence, it can be expected that 
with increasing HCB density beyond a certain value of 
the chemical potential a resonant state is occupied. This 
occurs when the expense in kinetic energy due to the loss 
of some degrees of freedom is compensated by a smaller 
interaction energy. Above the critical doping a HP can 
cluster with another on the same rung but then it tun- 
nels from one suitable HP to another HP thereby gaining 
additional kinetic energy. A set of parameter values for 



is given in Table [II. With these we calculated 



the HCB density per rung, (Ni) = Y'.j-i (-^i.?), for 3 up 
to 6 HCB on a 15 X 2 ladder shown in Fig. These 
results correspond to the hole density profiles shown in 
Fig. I and Fig. | for the 16 x 4 t-J ladder for 3 up to 6 
HP. As can be seen, the formation of FHC is reproduced 
qualitatively. The number of maxima in the density pro- 
files does not change, if a fourth HCB (HP in the t-J 
model) is added and for 6 HCB (HP) we observe three 
well separated maxima, each containing two HCB (HP). 
The density profiles for 5 HCB and HP look a bit differ- 
ent. Here the two additional HP cluster with the others 
but, due to the repulsive interaction between them, they 
are pushed towards the open ends of the sample in x- 
direction. We observe the same effect in the effective 
model, where the two outer maxima have higher weights 
than the maxima in the center. We tested the effective 




FIG. 13. Hard core boson (HCB) density per rung {Ni) 
obtained from the effective model (Il7|) for a (15 x 2) ladder 
with 3 up to 6 HCB. From below: 3, 4, 5, and 6 HCB, re- 
spectively. The plots are shifted by 0.2 with respect to each 
other. The interaction parameter values from Table HI have 
been used. The data reproduces qualitatively the hole cluster 
formation shown in Fig. ^ and Fig. ^ 



model also for longer systems. Figure 14 shows the rung 
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FIG. 14. Hole density per rung (n^) on a (24 x 4) t-J 
ladder with J = 0.35 1 for 5 up to 8 hole pairs doped into the 
ladder. From below: 5, 6, 7, and 8 hole pairs, respectively. 
The plots are shifted by 0.2 with respect to each other. 



density profiles for 5 up to 8 HP on a (24 x 4) t-J lad- 
der. The FHC start to form when going from 5 to 6 
HP doped into the ladder, which corresponds again to a 
critical doping ^ ~ i. By filling the ladder further with 
HP the number of maximum can also decrease. For 7 
HP we observe only four maxima. The area attached to 
each maxima suggests that two FHC form near the open 
ends by occupying seven rungs each, whereas the remain- 
ing three HP share the inner ten rungs, where one pair 
is loosely bound with the other two. Finally, for 8 HP 
we observe 4 FHC which occupy 6 rungs each. Also this 
"contraction" of HP can be reproduced with the effective 



model using the interaction parameters from Table HI. 
Figure [ij displays the corresponding rung density pro- 
files for 5 up to 8 HCB on a 23 x 2 ladder. They show 
the same evolution with HCB number as the hole density 
profiles in Fig. |lj with HP number. We cannot expect 
the effective model to reproduce the behavior of the t- 
J model in all details, but in view of the simplicity of 
the model, the results are satisfactory. After scanning a 
wide parameter range for v(rj;,ry) we came to the con- 
clusion that the formation of FHC can occur when the 
interaction leads to resonant states. 

Since we observe a qualitative change in the density 
profiles, we have computed the chemical potential ii{h) 
of the HP as a function of hole doping 



Ai(/i) = E{h+1) - E{h-1) 



(18) 



for both models in order to check whether it shows a sin- 
gularity at the critical doping. Here, E{h) denotes the 
ground state energy for a system with h holes or h/2 HCB 
in the t-J and effective model respectively. Figure |l^ 
displays the chemical potential for HCB in the bosonic 
model for a 23 x 2 ladder as a function of the hole dop- 
ing S in the corresponding t-J model. The hole doping 
is related to the HCB number Nt hy 6 = Nt/2L. The 
chemical potential shows a deviation from the quadratic 
form observed at low densities right at the position where 
the formation of the FHC sets in, i. e. when Nj, is in- 
creased from 5 to 6. The exact position of this kink can 
be tuned by choosing appropriate interaction parameter 
values v(rx,ry). The difference in the chemical potential 
between the critical and zero doping, A/i ~ 0.42 1, is to a 
good approximation given by the difference between the 
energies of the first resonant state and the ground state, 
AE = 0.43 1, obtained by solving the two particle prob- 
lem. The critical doping can be shifted to lower values by 
making the interaction longer ranged. Figure |l^ shows 
the chemical potential for hole pairs obtained from the 
corresponding 24 x 4 t-J ladder. The dashed dotted line 
shows a fit to a second order polynomial for S < Sc- As 
expected, the chemical potential /i(ft.) deviates from this 
polynomial for 6 > Sc- 

Finally, we note that while the effective model repro- 
duces the evolution of the system from HP to FHC nicely, 
it does not reproduce the density profiles quantitatively 
in the FHC density regime as can be seen, for example, 
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FIG. 15. Hard core boson (HCB) density per rung (A*';) 
obtained from the effective model dlTl) for a (23 x 2) ladder 
with 5 up to 8 HCB. From below: 5, 6, 7, and 8 HCB, re- 
spectively. The plots are shifted by 0.2 with respect to each 
other. The interaction parameter values from Table HI have 
been used. The data reproduces qualitatively the hole cluster 
formation shown in Fig. p^ . 
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FIG. 16. Chemical potential jx/t for hard core bosons 
(hole-pairs) obtained from the bosonic model ( [l7| ) for a 
(23 X 2) ladder using the interaction parameters from Ta- 
ble HI as function of the hole doping 5 in the corresponding 
t-J model. The dash-dotted line shows a fit to a second order 
polynomial in the region 5 < 5c- 
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FIG. 17. Chemical potential /i/t for hole-pairs as a func- 
tion of hole doping 5 obtained for a (24 x 4) ladder with 
J = 0.35 1. The dash-dotted line shows a fit to a second order 
polynomial in the region 5 < 5c- 

by comparing the upper curves in Fig. ^ and |l^. There- 
fore, further refinement of the effective model would be 
required to obtain the accuracy necessary to calculate 
parameters such as Kp from the effective model in this 
density regime. 



V. CONCLUSIONS 

In this work our previous study of the low energy prop- 
erties of the hole doped t-J ladders with two legs is ex- 
tended to the case of four legs. In both cases density 
matrix renormalization group (DMRG) results show that 
the holes bind in pairs at low densities and a finite spin 
gap is preserved. To analyze the DMRG results in more 
detail we introduce a hard core boson model on a single 
chain to describe the low energy degrees of freedom. The 
effective interactions between the hard core bosons which 
represent hole-pairs, are determined by fitting the den- 
sity profiles to those obtained by DMRG methods for the 
t-J model. This effective repulsive interaction is longer 
ranged for the four-leg ladder and this in turn reduces 
the value of the Luttinger liquid parameter Kp. As a 
consequence the region of predominantly superconduct- 
ing correlations is reduced in the wider ladder. Whereas 
the hole-pairs in the two-leg ladder simply repel each 
other, in the wider four-leg ladder a modification of the 
hole density profiles appears beyond a critical hole dop- 
ing which can be simply interpreted as the formation of 
four hole clusters. To reproduce this behavior in the hard 
core boson model, it is necessary to introduce an extra 
transverse degree of freedom by replacing the single chain 
with a two-leg ladder. It is also necessary to modify the 
interaction potential to incorporate a four hole cluster as 
a finite energy resonance. This is achieved by replacing 
the monotonically decreasing repulsive interaction on the 
single chain by one with a potential minimum at short 
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FIG. 18. Smoothing function /(d). The parameter d gives 
the distance to the nearest outermost rung of the system in 
units of the width m of the boundary. 



range. The longer range repulsive tail of the interaction 
potential causes the chemical potential to add a hole-pair 
to rise as the density is increased. When the chemical po- 
tential exceeds a critical value approximately equal to the 
resonance energy, the four hole clusters are formed lead- 
ing to a kink in the chemical potential as a function of 
hole density. 
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fid) = { 





1 

2 
1 



d=0 

d-i 1 

1 -I- tanh < < 1 

d>=l . 



(Al) 



and plotted in Fig.p^. Here d denotes the distance to the 
nearest outermost rung of the system in units of m. For 
the t-J model we use 

di,i' = dx^y,x',y' = min( ^ - 1, i - ^ )/m (A2) 

with L as the number of sites in the long direction of the 
t-J ladder and scale the energy parameters in such a way 
that tij/t = /(rfi.j) and Ji^j/J = /(dij), where t and J 
arc the bulk values. 

For the one dimensional chain with length N of the 
effective model we define 

, . ,R + R' ,^ R + R' . , / , ON 

dR^R, = mm( — 1,N — )/m (A3) 

dR ^ mm{R - 1,N - R)/m 

and scale the energy parameters in the way that t*j^/t* ~ 
/(rffl+i), ei?/e = /(di?), Wb.ii/wb = /(^i?), and 
Wint,fl,fl'/'^int = fidR^R/). Here again, t*, e, Wb, and Vi^t 
are the bulk values. 

To determine an appropriate value for the width m of 
the boundaries we have considered the curvature of the 
rung hole density profile of a (24 x 4) t-J ladder. For 
m < 4 the curvature is far from being a smooth function 
of the rung coordinate but for m = 4 it is. So we have 
used m = 4 for the computations with SBC. 



APPENDIX A: SMOOTH BOUNDARY 
CONDITIONS 

The SBC introduce smoothly decreasing energy pa- 
rameters into the Hamiltonian as the edges of the lattice 
are approached. [|l3|,0 The result of this operation is 
that, instead of having a sharp and rigid boundary, the 
boundary extends itself into the system and its exact size 
is not fully determinable. In this way we talk of the bulk 
of the system as the region where the energy parame- 
ters are constant, and of the boundary as the region over 
which the parameters are smoothly turned off. 

For the t-J ladders we use SBC only in the long direc- 
tion and denote the width of the boundary with m. By 
applying SBC to the Hamiltonians (|l|) and (^) wc replace 
the energy parameters t, J, <*, e, Vh, and Wint by the site 
dependent parameters iij, Jij, f^, er, Vh,R, and Vint,R,R' 
respectively. All these parameters are scaled according 
to the smoothing fimction / defined as 
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